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ABSTRACT 

The first investigations of tlie response of the microlensing magnification pattern 
(at an optical depth of order unity) to the mass function of the microlenses found 
that the resulting statistics depend only on the mean microlens mass (m). In par- 
ticular the mean microlensing caustic crossing rate was found to be proportional to 
a/ (m). We show that while this is true in the limit of mass functions with a narrow 
range of mass, in general the magnification pattern shows structure that reflects the 
contribution to the optical depth of microlenses with different masses. We present a 
better approximation, relating the microlens mass function to light-curve statistics. 
We show that the variability statistics of quasar microlensing light-curves can (in prin- 
ciple) be inverted to obtain the mass function of the microlenses in the mass range 
over which the mass density remains comparable, ie. p(m)dm ~ Cm^^. A prelimi- 
nary analysis of the structure function for Q2237-I-0305 suggests that there is not a 
significant contribution to the optical depth from very low mass objects {\Q~^Mq). 
However observations of multiple microlensed quasars for a period of ~ 20 years may 
in the future yield a detailed p(rn)dm. In the mass range where the number density 
is comparable, ie. p{m)dm « const., the distribution of flux factors could be inverted 
to find the microlens mass function. This may be used as a probe of the abundance of 
planets with orbital radii > lOOAU. 



Key words: gravitational lensing ■ 
Planets: masses. 
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1 INTRODUCTION 

Following the findings of Witt, Kayser & Refsdal (1993) and 
Lewis & Irwin (1995, 1996) the mean microlens mass has 
been assumed to be the only important mass function pa- 
rameter with respect to quasar microlensing statistics. The 
primary implication of this assumption for the determina- 
tion of the microlens mass function is that observations of 
quasar microlensing could yield a value for the microlens 
mass independently from the detail of the mass function. 
However the converse, that quasar microlensing provides no 
information about the detail of the microlens mass func- 
tion (Witt, Kayser & Refsdal 1993) is contrary to the initial 
hopes that such information could be extracted from long 
term monitoring of objects like Q2237-f0305 (e.g. Wambs- 
ganss, Paczynski & Schneider 1990). 

Lewis & Irwin (1995) showed that the point source mag- 
nification distribution is independent of the mass function. 
However for a given set of macro-parameters they also show 



that the magnification distribution is affected by a compo- 
nent of optical depth in smooth matter (kc). On the other 
hand it is not surprising that the magnification distribution 
depends on the optical depth (k) and the shear (7). The de- 
pendence on ACc is therefore readily understood in light of the 
parameter transformation of Paczynski (1986) which relates 
models containing continuous matter to those that include 
only point masses but are described by different macro pa- 
rameters. In combination with a macrolens model, the mag- 
nification distribution provides a means to probe the fraction 
of continuous matter (Lewis & Irwin 1995). Of course the 
effect of a finite source size, which narrows the width of the 
magnification distribution, also needs to be considered. 

Refsdal & Stabell (1993) noted qualitatively that if both 
small and large masses contribute significantly to the optical 
depth, one expects that, for a source size lying between the 
respective microlens Einstein Radii (ER), small fiuctuations 
will be superimposed on the large scale variations caused by 
microlensing due to the larger masses. This statement must 
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be true for the following reason. If the example is taken to its 
extreme, by decreasing the small mass and at the same time 
keeping the optical depth due to the small mass population 
constant (ie. by increasing their number density), then the 
model acts like one that includes point masses plus a con- 
tinuous sheet. The mean mass is zero, but the microlensing 
rate certainly is not. Cl early microlensing variability is not 
determined solely by yj' (m). 

In this paper we discuss the effect of a mass function 
on various microlensing statistics, and show that in principle 
the microlens mass function can be retrieved from long-term 
monitoring of microlensed quasars. Sec. ^ presents magni- 
fication patterns for different mass functions, and gives a 
qualitative discussion of their features. In Sec. |^ structure 
functions and flux factor distributions are shown. Approx- 
imate expressions are presented which relate the microlens 
mass function to the microlensing light-curve statistics. Fi- 
nally in Sees. ^ and ^ we present a preliminary analysis of 
the structure function for Q2237+0305, and discuss poten- 
tial applications. 



2 MAGNIFICATION PATTERNS 

In this section we present magnification patterns produced 
by microlens populations with different mass functions hav- 
ing upper (7712) and lower (mi) limits that produce mass 

m2 



^ = 100 and ^ = 1000. Four 
classes of mass function have been chosen for illustration: 



ranges oi^ ^1,^ = 10, 



i) Psingieim)dm = S{m) 

... , S{mi)+5{m2) 

II) PNdense{m)dm = 



*) PSalpeter{in)dm (X m 



nil < m < 1712 



otherwise 



7722 

iv) podepth(ni)dm oc S(mi) + 5(m2) 

mi 

Psingieim) describcs a population of identical point masses. 
PNtienseim) and podepth{fn) describe mass distributions with 
equal number densities and optical depths respectively in 
the two populations of masses mi and m2. psaipeter{m) is a 
Salpeter mass function. 

The ray-tracing method (e.g. Kayser, Refsdal & Sta- 
bell 1986; Schneider & Weiss 1987; Wambsganss, Paczyn- 
ski & Katz 1989) was used to calculate the magnification 
patterns. The number of stars used in the models was calcu- 
lated through the method described in Lewis & Irwin (1995) 
and Wyithe & Webster (1999). Fig. |l| shows 4 columns of 
magnification patterns, each 20 ERia/q (Einstein radius for 
IM0) on a side. The columns show magnification patterns 
corresponding to mass functions psingie{m), pNdense{m), 
psaipeter{m) and podepthim), each with an optical depth of 
K = 0.4 and no applied shear. The upper 4 rows show mag- 
nification maps produced by microlens populations with de- 
creasing (m) and increasing mass range (^). The mean is 
constant in each row and was determined so that m2 = IA/q 
in column 4 {podepth{m)): 



(1) 
(2) 

(3) 
(4) 



(m) 



(5) 



The bottom row shows the limiting case as mi — > 0. The 
number of stars and the mean magnification respectively for 
these magnification maps are shown in Tab. ^. 

Columns 2—4 of Fig. ^ demonstrate that the introduc- 
tion of different scales into the structure of the magnification 
map results from the inclusion of comparable mass-density 
from microlenses at opposite ends of the mass range. Fig hi 
shows several additional trends. The magnification maps in 
column 1 {psingie{m)) show the well documented increase 
in uniformity as m is lowered (e.g. Wambsganss, Paczyn- 
ski & Katz 1989). The magnification pattern approaches a 
constant for a small but finite source. This is also true for 
the other mass functions discussed. However, the magnifi- 
cation maps corresponding to PNdense{m), Psaipeter{m) and 
POdepthijn) in columns 2, 3 and 4 demonstrate a slowing of 
this trend if mass density is distributed over a finite mass 
range. For podepthim) where half the mass density is in IMq 
stars and the other half in 10~^Mq microlenses, there is no 
increase in uniformity over the single mass case, even though 
the area in square ER (of the mean mass) is increased by 



a factor of 500. Rows 3 and 4 (^ 



100, 1000) show a 



decrease in uniformity of the magnification map were con- 
tributions to the mass density of microlenses are more evenly 
spread across the mass range. Only a small section of mag- 
nification map is presented (in Fig. |^) for the smaller mass 
range examples. Magnification maps corresponding to the 

rper 4 rows of Fig. hi have therefore also been plotted in Fig. 
this time with a mean microlens mass of (m) = O.2A/0 in 
all cases (the simulation details are summarised in Tab. 
Fig. |2| shows the trends already mentioned and in addition 
demonstrates the more significant loss of uniformity where 
the mass range is greater. However, note that the magnifica- 
tion map for podepthim) clearly shows larger scale structure 
than the magnification map for psingieijn) even when the 
mass range is only 10 in the former case. 

Podepthim) may not be physically realistic, however 
the resulting magnification pattern demonstrates the major 
point of this paper. The magnification maps produced by mi- 
crolens mass functions of the form podepthim) appear to be 
averages of the magnification maps produced by psingieim2) 
and psingieim\). This is demonstrated by Fig. ^, and sug- 
gests that microlenses contribute to the magnification pat- 
tern in proportion to their mass density. We explore this 
statement in the next section. 



3 LIGHT-CURVE STATISTICS 

We investigate the relationship between the mass function 
and the magnification pattern using two light-curve statis- 
tics. The first is the structure function S(Aj/), defined by: 



S(Ay) = {lA/(y)-A//(y + A2/)|), 



(6) 



where y is the source position (or equivalently an epoch) , and 
Miy) is the corresponding image magnitude. The second 
statistic is the distribution of fiux factors piK)dK where K 
is defined by the near caustic approximation of Chang & 
Refsdal (1979): 



K 



(7) 
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Figure 1. Magnification patterns for different mass functions. Each row has a different mean microlens mass (shown in column 1), and 
each of the columns 1-4 a different mass spectrum (mass functions p Singles PNdensey PSalpeter ^i^d Podepth respectively). The upper and 
lower mass limits are shown by mi and m2 in each case. The bottom row shows the limiting case where — > 0. 



In Eqn. |7|, is the point source magnification, /lo is the mag- 
nification of non-critical images, Aycaust is the distance of 
the source from the caustic, and O is the Heaviside step func- 
tion. We calculate K using the expression of Witt (1990). 
Using S and p{K)dK, the statement at the end of Sec. H, 



relating the microlens mass function to the magnification 
pattern can be written as: 

S{Ay) = ^ J mp{m)SiM,, (^^^ dm (8) 
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Figure 2. Magnification patterns for different mass functions, with a mean microlens mass of (m) = 0.2Mq in all cases. The columns 
1-4 show results from different mass spectra (mass functions psingle^ PNdeiiae^ PSalpeter a-nd POdepth respectively). The upper and lower 
mass limits are shown by mi and m2 in each case. 



Table 1. Values for the number of stars used and the average magnification obtained from the simulations presented in Fig. y. The 
theoretical mean is /lav = 2.78. 



(m) (Mq) 


m2 
mi 


P Single 


Mass function 

PN dense PSalpeter 


POdepth 


1.0 


1 


592 (3.24) 


592 (2.34) 


592 (2.93) 


592 (2.99) 


0.18 


10 


1970 (2.91) 


2222 (2.44) 


2165 (2.56) 


2637 (2.82) 


0.02 


100 


13308 (2.74) 


14193 (2.80) 


15643 (2.81) 


23576 (2.78) 


0.002 


1000 


117599 (2.78) 


120284 (2.79) 


136754 (2.79) 


233018 (2.85) 


NA 


NA 


NA (2.79) 


NA (2.79) 


NA (2.79) 


232 (2.24) 
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Figure 3. Magnification patterns for single mass microlens populations (psingle) of m = IMq (left) and m = O.OIM0 (centre), and 
for a model with podepth having = 100 and m2 = 1 (right). The right-hand panel is approximately an average of the left-hand and 
centre panels. 



Table 2. Values for the number of stars used and the average magnification obtained from the simulations presented in Fig. Q The 
theoretical mean is fiav = 2.78. 



Mass function 



(m) (Mq) 


m2 
mi 


P Single 


VN dense 


PSalpeter 


POdepth 


0.2 


1 


117599 (2.77) 


117599 (2.77) 


117599 (2.77) 


117599 (2.77) 


0.2 


10 


117599 (2.77) 


119495 (2.78) 


119084 (2.77) 


122428 (2.72) 


0.2 


100 


117599 (2.77) 


120202 (2.79) 


124354 (2.79) 


145246 (2.82) 


0.2 


1000 


117599 (2.77) 


120284 (2.79) 


136754 (2.85) 


233018 (2.79) 



Km 



dm. 



(9) 



and 

p(K) = J ^/mp{m)piMQ 

Here Simq and PiMq are the structure function and distri- 
bution of flux factors for a population of IMq microlenses. 
p{m)dm is the microlens mass function and Ns and Nk are 
normalising constants such that: 

Ns = J mp{m)dm (10) 
and 

Nk = J y/mp{m)dm. 

Eqn. ^ contains an extra factor of since the length of 
caustic per cm'^ increases with -j^- For comparison, the as- 
sumption that light-curve statistics only on the mean may 
be written as: 



(11) 



5'(Aj/) = Su 



and 



Ay 



(12) 



(13) 



respectively. 

In the limits of small mass ranges (mi m^) or where 
the dominant fraction of optical depth is present in mi- 
crolenses of a single mass mo (ie. mp{m) moS{mo)), Eqns. 



Table 3. Values for the number of stars used, and the average 
magnification obtained from the light-curve simulations. The the- 
oretical mean is /lav = 2.78. 







Mass function 




mi 


PN dense 


PS alpeter 


POdepth 


1 


1644 (2.73) 


1644 (2.73) 


1644 (2.73) 


10 


1833 (2.72) 


1791 (2.64) 


2153 (2.70) 


100 


1907 (2.67) 


2379 (2.70) 


5622 (2.57) 



and ^ reduce to Eqns |l2| and In addition Eqns. ^ and 
are correct in the limit of low optical depth. If Eqns. ^ 
and ^ describe the relationship between the magnification 
map and the microlens mass function, then these and simi- 
lar expressions offer the means to measure the mass function 
from microlensing light-curve statistics. In the next section 
we explore this idea using model microlensing light-curves. 



3.1 the microlensing models 

We have investigated structure functions and distributions 
of flux factors produced by models with m = IMq for mass 
function psingie and by models with (m) = IMq for mass 
functions Pn dense, PSalpeter aud po depth with ^12.^10 and 
^ = 100. For each model, 500 point source light-curves 
were produced (one per random starfield), each with a length 
of 60 ERimq using the contouring method of Lewis et al. 
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Magnitude change Magnitude change Magnitude change 

Figure 4. Magnification distributions for models with = 1 (solid lines), = 10 (dashed lines) and = 100 (dotted lines). 
Left: Distributions for mass function Pn denser Centre: Distributions for mass function Psalpetery Right: Distributions for mass function 

POdepth- 



(1993) and Witt (1993). The number of stars in each model 
was determined from the description of Lewis & Irwin (1995) 
(using the formalism of Katz, Balbus & Paczynski (1986)). 
Table ^ shows the number of stars and the mean magnifica- 
tion for each model together with the theoretical magnifica- 
tion. The resulting magnification distributions are shown in 
Fig. 0. Note that we confirm the finding of Lewis & Irwin 
(1995) that the magnification distribution is independent of 
mass function. 

3.2 the structure function 

Eqn. 1^ suggests that the structure function has a shape that 
reflects the microlens mass distribution. Small masses pro- 
duce more rapid variability and therefore a faster rise in 
the structure function at small Ay. Larger masses cause 
the asymptotic behaviour of the structure function at large 
Ay to slow. Fig. ^ shows structure functions corresponding 
to microlens mass functions of the form pNdense, psaipeter, 
POdepth, with ^ = 100. These demonstrate the effect of a 
varying mass function on the structure function of a mi- 
crolensed light-curve, as well as the applicability of Eqn 
^ In each case 4 curves are shown in the left-hand figure. 
The solid light and solid dark lines correspond to structure 
functions computed directly from models and from Eqn. ^ 
respectively. For comparison, the structure function corre- 
sponding to the single mass microlensing model (psinaie) is 
shown by the dot-dashed line. Fig. ^ shows that Eqn. ^ pro- 
vides a significantly better approximation to the microlens- 
ing statistics than Eqn. |l^. To quantify the success of the 
approximation we have fitted for the three free parameters 
of p(m) in Eqn. ^ ie. (m), ^ and p{m\) or the index a (for 
fits corresponding input mass functions Pn dense and po depth 
or psaipeter respectively). The input mass weighted number 
densities (those used to compute the model light-curves) 
are shown as the light histograms in the right hand pan- 
els of Figs. ^ In each case, the fitted mp{m) is overlayed 
(dark histogram). The fitting procedure, using Eqn. |^ is far 
more successful for mass functions which provide significant 
fractions of optical depth over a large mass range. This is 



because the dimensions of the magnification pattern only 
scale with y'm (so that very different masses are needed to 
effect the scale of the magnification map), and also because 
(from Eqn. ^ the mass function contributes to the differ- 
ent scales {\/rn) of the caustic network in proportion to the 
optical depth associated with mass m (ie. K(m)dm). The in- 
put and fitted parameters for pirn) with ^ = 10, 100 are 
summarised in Tab. ^ 

Fig. 1^ shows excellent agreement between the directly 
calculated, predicted and fitted structure functions in the 
case oi pm dense- However the mass function corresponding to 
the best fit is not in agreement with the input mass function. 
The reason is that the mass weighted number density for 
PNdense IS approximately that of mass function psingie with 
m = m2. Hence S{Ay) —> Simq [Ay / y/m,2). The contribu- 
tion of the small mass to the shape of the structure function 
is limited to small Ay, and S less than S{Ay — > cxd) x ^i^iii. 
The fit has a mean mass near m2 with a range smaller than 

m2 
mi 

Fig. 1^ shows reasonable agreement between the di- 
rectly calculated, predicted and fitted structure functions 
for podepth- In addition, unlike the case of PNdense, the 
fitted mass function is also in reasonable agreement with 
the input mass function. The contributions to the structure 
function from the two mass populations are distinct. Below 
Ay ~ 0.5ER, S is dominated by the small mass popula- 
tion which causes rapid variability. Above Ay ~ 0.5ER, S 
is dominated by the larger masses responsible for the longer 
term variability. There is a direct correspondence between 
two sections of S, and the two scale behaviour evident in the 
corresponding magnification map (Figs. ^ and ^). Because 
of the break at Ay ~ 0.5, S does not have a shape that is 
similar to a scaled version of Sm- The fit for p{m) is therefore 
much more successful. 

Finally, Fig. |^ also shows reasonable agreement between 
the directly calculated, predicted and fitted structure func- 
tions for psaipeter- The shape of S'(Ay) differs from that 
of S'lAfg (Ay/-^/{m)) because the smaller masses produce a 
more rapid rise to the structure function, while the larger 
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Figure 5. Left: The microlensing structure function. The hght, dark and dashed lines show the calculated and predicted structure 
functions and the best fit respectively. The curve corresponding to Psingle shown for comparison (dot-dashed line). Right: The input 
mass function (light line) and the fitted mass function (dark line). Results are shown corresponding to mass functions dense (top), 
PSalpeter (centre), and podepth (bottom). The bars in the upper and lower right-hand panels should be interpreted as having the same 
linear width. 



masses cause the asymptotic behaviour to slow. The input 
index a of psaipeter is approximately recovered. 

Our modelling is restricted to mass ranges of ^ = 100 
or smaller, and the effect of the mass function is subtle in 
all but a couple of cases. This is because of the -^m depen- 
dencies of microlensing statistics. Mass functions with much 



larger mass ranges would show more significant effects. How- 
ever the results of this section demonstrate that Eqn. ^ is a 
quantitative improvement over Eqn. |l^, and that it is qual- 
itatively correct in many cases. 
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Table 4. Summary of input and fitted parameters for p{m) obtained from Eqn. |8| as well as the for the fit. 



true parameters 



fitted parameters 



Mass Function 


m2 
mi 


(m> 


/ / " 


m2 
mi 


(m) 


//" 


(128 points) 


PN dense 


10 


1.0 


0.5 


5.15 


1.13 


0.17 


7.2 


X 


10" 


-4 


PNdense 


100 


1.0 


0.5 


2.05 


1.81 


0.37 


1.2 


X 


10^ 


-3 


PSalpeter 


10 


1.0 


-2.35 


576 


1.00 


-3.06 


7.0 


X 


10" 


-4 


PSalpeter 


100 


1.0 


-2.35 


186 


0.76 


-2.03 


1.0 


X 


10" 


-3 


POdepth 


10 


1.0 


0.091 


29.5 


0.48 


0.03 


8.4 


X 


10" 


-4 


POdepth 


100 


1.0 


0.010 


281 


0.41 


0.003 


2.5 


X 


10" 


-3 




^.01 





0.1 



1 

K 



10 



100 



Figure 6. Distributions of flux factors for microlensing by a two mass population (pNdenae)- The light, dark and dashed lines show 
the calculated and predicted distributions and the best fit respectively. The curve corresponding to (p Single) shown for comparison 
(dot-dashed line). Right: The mass function for PNdense. (light line) and the fitted mass function (dark line). 



3.3 the distribution of flux factors 

The previous section demonstrated the dependence of mi- 
crolensing variability statistics on the microlens mass func- 
tion. The size of a light-curve peak due to a caustic crossing 
is dependent on the flux factor K (a function of microlens 
mass). In addition to the dependence of the variability, it 
is therefore reasonable to suppose that the distribution of 
peak heights should be related to the mass weighted num- 
ber density, with an additional weighting of to account 
for the higher caustic density (Eqn. |^. 

The same set of simulations and fits discussed in Sec. 
^ were made for p{K)dK, and the results are summarised 
in Tab. ^. The extra factor of in Eqn. |^, and the de- 
pendence on rather than y'm means that p{K)dK and 
S{Ay) carry different information. In particular, for the case 

of equal optical depths in two disparate mass populations 

_i 

p(K)dk pm{Km-^ as milmx becomes large. As a con- 
sequence, in contrast to flts using S, the fits for psaipeter and 
POdepth are poor. However, the factor of boosts the con- 
tribution of the small masses in pNdense , results for which are 
shown in Fig. || (for ^ = 100). The calculated, predicted 
and fitted functions p{K)dK show excellent agreement, and 
the fitted mass function is also good. The reason for the 
more accurate recovery of the mass function is the bump in 
the distribution at small K corresponding to caustics related 
to the small masses. 

Fig. R shows two predicted flux factor distributions, cal- 



culated using Eqn.|9|. The left hand figure shows the distribu- 
tion produced by pNdense, with m2 — IMq and ^ — 1000 
(row 4 of column 2 of Figs. ^ and ^). This distribution 
shows a "foot" which is the signature of the small (~Jupiter 
mass) microlenses. The right hand figure shows the distri- 
bution produced by a two population mass function with 
p(mi) = 10p(m2), m2 = IMq and ^ = 100000. The very 
low (~Earth) mass objects contribute a spike feature to the 
p{K) near K ^ 0. 



4 PRELIMINARY APPLICATION TO 
Q2237+0305 

Lewis et al. (1996) found that the structure function for 
Q2237+0305 (Huchra et al. 1985) suggests a mean mass of 
O.li\f0 ^ (m) ^ 1OM0. This finding remains true in light of 
additional light-curve data. There are currently ~ 190 light- 
curve points for each of the 4 images of Q2237+0305, span- 
ning a baseline of ~ 5000 days (Schneider et al. 1988; Kent & 
Falco 1989; Irwin et al. 1989; Corrigan et al. 1991; 0stensen 
et al. 1996; Wozniak et al. 2000a,b). In the past three sea- 
sons, the OGLE collaboration (Wozniak et al. 2000a,b) has 
obtained monitoring data at a high-rate and with good 
photometric accuracy. For the first time this data probes 
short duration fiuctuations. Fig. ^ presents the binned struc- 
ture function, averaged over the 4 images of Q2237-(-0305. 
Quoted errors in the data have been added in quadrature for 
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Figure 7. Predictions of distributions of flux factors (Eqn. |sj) for microlensing by a two mass population. Left: The distribution for 
mass function pf^ dense with = 1000. Right: The distribution for a two mass population with = 100000, and p{mi) = 10p(m2). 
The curve corresponding to psingle. is shown for comparison (dot-dashed line). 



Table 5. Summary of input and fitted parameters for p{m) obtained from Eqn. bl as well as the for the fit. 



true parameters fitted parameters 



Mass Function 


mi 
m\ 


(m) 


//" 


m2 
Ta\ 


(m) 


/ / " 


(128 points) 


// 


10 


1.0 


0.5 


6.30 


0.89 


0.11 


4.6 X 10-2 


// 


100 


1.0 


0.5 


66.90 


0.59 


0.01 


0.1 


/// 


10 


1.0 


-2.35 


5.16 


0.92 


-0.75 


5.6 X 10-2 


/// 


100 


1.0 


-2.35 


783 


0.99 


-3.37 


5.5 X 10-2 


IV 


10 


1.0 


0.091 


5.15 


0.906 


0.17 


5.5 X 10-2 


IV 


100 


1.0 


0.010 


9.53 


0.863 


0.10 


9.6 X 10-2 



each pair, and subtracted from the square of the magnitude 
difference before binning. The error bars for the structure 
function were estimated as the variance in each bin (Af) di- 
vided by the squareroot of the the monitoring period over 
or the number of points, whichever is less. 

Superimposed on these plots are the model structure 
functions computed for a source size of 0.13 ERo.ia/q and 
a transverse velocity of 400 fern see" ^. The solid dark line 
shows the structure function for podnpth with m2 = O.IMq 
and ^ = 100. This was computed using a magnification 
map (k = 0.4, 7 = 0) of side-length 20 ERimq and com- 
posed of 500 X 500 pixels. The dashed line shows the corre- 
sponding structure function computed to smaller At using 
Eqn. 1^ and higher resolution single mass (m) microlensing 
magnification maps (side- length IO^^ERimq and com- 
posed of 1000,/^ X 1000,/^ pixels). The light solid line 
shows the structure function for psingie with m — O.IMq 
computed from the 1000 x 1000 pixel magnification maps. 

In Wyithe, Webster & Turner (2000) it was demon- 
strated that the source size must be smaller than ~0.2 ERmj 
from the large scale variability of recent seasons. The as- 
sumption of a source smaller than 0.13 ERmi will result 
in more rapid initial rise of the structure function for the 
model with podepth- Note that for this model, the source is 
~ lERmj^ , and that the asymptotic value is slightly lower, 
which corresponds to a narrower magnification distribution. 
Both model structure functions approach their asymptotic 
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Figure 8. Binned structure function averaged over the 4 images 
of Q2237-I-0305. Superimposed are model structure functions for 
both a single mass model {p single) with m = O.IMq (light line) 
and a podepth model with m2 = O.IMq and = IQO (dark 
line). The dashed line is an approximate structure function for 
the Podepth model computed using Eqn. ^. The assumed effective 
transverse velocity was 400fcmsec-^, and the source size 0.13 
ERmo • 
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value at around 5000—10000 days, reflecting the microlens- 
ing time-scale of the 0.1 Mq microlenses. 5000 days is ap- 
proximately the length of the monitoring period, and is 
therefore the separation where the observed structure func- 
tion is least determined. However, the structure function is 
better determined at shorter separations. In particular if the 
microlens population in Q22374-0305 can be described by 
POdepth with 1712 = O.IM0 and ^ — 100 (dark solid line), 
then the effective sampling length for the structure function 
at t ^ 500 days is increased by a factor of 10, so that it be- 
comes representative on these time-scales. However qualita- 
tive comparison between the models and the structure func- 
tions for Q2237-)-0305 does not show evidence for a contri- 
bution to the optical depth from very low mass microlenses. 
A proper treatment of this result will require models with 
appropriate values for k and 7, as well a monte-carlo analy- 
sis of structure functions modelled using sampling rates and 
simulated errors corresponding to those of the observations. 
Such an analysis will be presented in a future paper. 



5 DISCUSSION 

The microlensing magnification pattern has a structure that 
reflects the mass spectrum of microlenses. In particular, vari- 
ability statistics for a given, mass function are are well ap- 
proximated by the average of variability statistics for single 
microlens mass distributions, weighted by the correspond- 
ing microlens mass and number density. Put another way, 
microlenses of different masses contribute to the variability 
statistics in proportion to n{m)dm. In addition. High Mag- 
nification Event (HME) statistics for a particular mass func- 
tion are well approximated by the average of HME statistics 
for single microlens mass distributions, weighted by the mi- 
crolens mass (m), the number density, and where the 
last factor is included to account for smaller masses produc- 
ing more caustics per unit area. 

If there is a similar contribution to the optical depth 
from masses at opposite ends of the mass range, the variabil- 
ity statistics are non-degenerate with the single mass case, 
and the mass function can be inferred by inverting an inte- 
gral (Eqn. We have carried out simple parametric fits for 
the mass function using models with known distributions to 
test this idea. In practice however, more sophisticated, non 
parametric inversions could be used to obtain the mass func- 
tion from variability data. On the other hand, since Eqn. ^is 
an approximation, an inversion will most readily be achieved 
by directly comparing computed model statistics with data. 

The distribution of flux factors describes the size of 
HME peaks. The distribution shape is sensitive to a different 
set of mass functions than the structure function. In particu- 
lar, since smaller masses contribute more caustic length per 
unit area, the statistics of HMEs provide information on the 
mass function where the number density of objects at dif- 
ferent ends of the mass range, rather than the optical depth 
is comparable. Brown dwarfs or planetary mass objects con- 
tribute to the flux factor distribution as an excess of small 
values. Our models assume random distributions of objects, 
and it is not clear what effect the spatial correlation of the 
large and small masses will have on the magnification map. 
However, in projection, and at optical depths of k ~ 0.5, 
planets with orbital radii > lOOAU are not spatially corre- 



lated with stars. Microlensing of quasars by galaxies at mod- 
erate redshift may therefore shed light on the abundance of 
planets in a regime not accessible by other methods. 

Our investigation has concentrated on what informa- 
tion about the distribution of microlens masses is mani- 
fest in magnification maps, and in the corresponding light- 
curve statistics. However, in practise obtaining sufficient 
quality data will be difficult. Currently, only one object, 
Q2237+0305, is known to exhibit significant quasar mi- 
crolensing. However, the Sloan Digital Sky Survey will ob- 
tain spectra of 10^ galaxies (York et al. 2000), and it is 
estimated that a few 10s of these may have very close align- 
ment with a background quasar, as well as a small redshift, 
and will therefore be efficient microlenses like Q2237-I-0305 
(Fukugita & Turner, Private Communication). With a large 
number of monitored objects, statistics of microlensed light- 
curves will become much better known. 

We have computed statistics for a point source and 
given transverse velocity, while in practice both the source 
size and transverse velocity will be unknown. However tech- 
niques have been developed to estimate these quantities from 
monitoring data (Wyithe, Webster & Turner 1999,2000; 
Wyithe, Webster, Turner & Mortlock 2000). Long light- 
curves extending over 10s of Einstein Radii for a IMq 
star will not be available for a century or more. However 
the structure function effectively levels out after ~ 2 mi- 
crolens (mean) Einstein Radii. Therefore, for typical galac- 
tic proper motions, ~ 20 years of monitoring for multiple 
objects should yield an observable signal. 

The detailed shape of the mass distribution is contained 
in the combination of a structure function and distribution 
of flux factors of sufficient quality. However the structure 
function can be used to place limits on the fraction of opti- 
cal depth in small objects from the initial slope. Moreover, 
the time scale at which the structure function levels out 
yields the upper limit. The more moderate goal of probing 
the mass range or the existence of two distinct mass pop- 
ulations should therefore be achievable. Similarly, limits on 
the existence of small mass objects will come from the lack 
of observation of small light-curve peaks. 
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